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Abstract: We study the dynamics of a strongly-coupled quantum field theory in a cosmo¬ 
logical spacetime using the holographic AdS/CFT correspondence. Specifically we consider a 
confining gauge theory in an expanding FRW universe and track the evolution of the stress- 
energy tensor during a period of expansion, varying the initial temperature as well as the rate 
and amplitude of the expansion. At strong coupling, particle production is inseparable from 
entropy production. Consequently, we find significant qualitative differences from the weak 
coupling results: at strong coupling the system rapidly loses memory of its initial state as 
the amplitude is increased. Furthermore, in the regime where the Hubble parameter is much 
smaller than the initial temperature, the dynamics is well-modelled as a plasma evolving 
hydrodynamically. 
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1 Introduction 

The dynamics of quantum fields in curved spacetime is of immense physical importance in our 
universe. Quantum fluctuations and particle production provide the seeds for the macroscopic 
structures we observe in the present epoch. While the standard inflationary paradigm works 
extremely well to explain the observations to date, it is based on the dynamics of free fields. 
On the other hand, it is conceivable that strongly-coupled field theory dynamics played a 
significant role during some phase of the very early universe. This motivates the study of 
strongly-coupled QFTs in cosmological spacetimes. 

Such an investigation would be impractical using direct field theory methods, which 
at best can quantitatively answer only very simple questions when strong interactions are 
present. However, if the field theory in question is holographic - i.e., is dual via the AdS/CFT 
correspondence to a gravitational theory in an asymptotically Anti-de-Sitter spacetime - 
we can translate questions about strongly-coupled cosmological physics into much simpler 
questions in classical gravity. This approach provides a powerful quantitative insight into 
dynamics of strongly-coupled QFTs in curved spacetime; cf., [1] for a review. 

In this paper, we focus on a classic question in the study of quantum fields on curved 
spacetime. Starting from a thermal equilibrium state in Minkowski space, we investigate the 
final state of the field theory after a period of homogeneous and isotropic expansion (see [2] for 
other scenarios). In general the evolution entails time-dependent couplings in the field theory 
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which lead to energy and entropy production. We are interested in ascertaining if there are 
qualitative differences in this production when compared with results in weakly-coupled field 
theory, where it can be understood as particle production in a time-dependent background. 
Roughly speaking, at strong coupling particle production is expected to be accompanied by 
rapid local thermalization, a feature that is postponed at weak coupling to the end of the 
expansion. We will see that this lack of separation of time scales in the strongly coupled 
theory leads to qualitatively different results from the weakly coupled context. 

The simplest field theories to study holographically are certain strongly-coupled con¬ 
formal field theories (CFTs) with many degrees of freedom. However, a CFT on an FRW 
cosmological spacetime does not have particle/energy production; a FRW spacetime is related 
by a conformal transformation to a static universe . 1 We therefore consider a non-conformal 
QFT 3 + i obtained by compactifying a CFT 4 + i on a circle. Studying this QFT 344 on an FRW 
spacetime is equivalent to studying the CFT 4 + i on a spacetime B = FRW x S 1 . 

Using the holographic duality, this translates to a classical gravity question of finding 
asymptotically AdSe spacetimes whose boundary geometry is B. By solving the relevant 
partial differential equations numerically, we can find solutions corresponding to a given initial 
temperature and FRW scale factor a(t). The physical data of the final state (thermodynamic 
and otherwise) after the expansion can be read off from the solution. 

Using these methods, we study how the amount of energy/entropy production are related 
to the rate and duration of cosmological expansion. Comparing with free held theory, we find 
interesting differences. As the amplitude is increased, the amount of energy produced during 
the expansion quickly becomes independent of the initial temperature, a result not seen in the 
free held limit. Further, at strong coupling when the maximum Hubble parameter is much 
smaller than the scale set by the temperature, the results for the temperature change show 
excellent agreement with analytic results from a hydrodynamic model, providing a strong 
check of the numerical results. 

2 Generalities 

We consider a holographic CFT 4 + i on the spacetime 

ds 2 = —dt 2 + dx 2 + a 2 (t) dx 2 , (2.1) 

where x c parameterizes a S 1 of size £ c , and anti-periodic boundary conditions are imposed 
for fermions. The scale factor a(t) evolves from a = 1 in the limit t = —00 to a = a at t = 00. 
For explicit calculations, we choose 

, , a + 1 a — 1 

a(i) = —- 1 -— tanh (vt). ( 2 . 2 ) 

1 Of course in making this assertion we are eliding over the fact that for even spacetime dimensions, CFTs 
have conformal anomalies. The conformal anomaly leads to a coupling independent Casimir type energy which 
is simply given by the local scale factor. It does not thus capture the physical effect of particle production we 

are after. Note however that it does have physical consequences for discussion of inflationary physics, cf., [3]. 
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This resembles a quench occurring around t = 0 and lasting for a duration v _1 . The compact 
S 1 provides an energy scale £f 1 to the effective 3+1 dimensional theory. At fixed scale 
factor (i.e., on Minkowski spacetime), this theory is in a confined phase for temperatures 
T < T c ~ £~ 1 * and in a deconfined phase for higher temperatures. We will focus, for simplicity, 
on the case where the field theory remains in its deconfined phase during the entire evolution. 
Thus we have some large number c e ff of deconfined degrees of freedom potentially participating 
in the process of particle production in the expanding background. 

We will explore how the final energy, entropy, and temperature depend on the initial 
temperature To, the amplitude a, and the rate v of the expansion. 

In the deconfined phase at large c e fr, local 4+1 dimensional quantities do not depend on 
the scale £ c , a property known as “large N volume independence” [4]. Holographically this 
follows from the fact that the finite £ c solutions are obtained from the £ c = oo solutions by a 
trivial identification x c ~ x c + £ c . Using this and scale invariance we we can parameterize the 
final 3+1 dimensional energy density in terms of the dimensionless quantities a and b = v/T 
as 

e(a,v,T 0 ) = C£ c v 5 f e (a, o) . (2.3) 

where the normalization constant C is related to the equilibrium energy density by eo(To) = 
C£ c Tq. The function / e (a, d) encodes the non-trivial dynamical information. 

In the limit v —> 0 of adiabatic expansion, the final energy density can be determined from 
the initial energy density using entropy conservation. The final state is related to the initial 
state by simple dilution and red shift effects. To isolate the effects of particle creation/energy 
production, we can compare the final state data relative to that attained simply during the 
adiabatic expansion for the same starting values. We thus define 

T e (a,t>) = / e (a,t>) - / e (a,D -+ 0). (2.4) 

Alternatively, we can consider the change in temperature, relative to the adiabatic result. 
Using conformal invariance, we find 

Tj = Ty(a, d — y 0 ) + v Gt(o, d) , (2-5) 

where the first term gives the adiabatic result, and Gt(c i, t>) contains the non-trivial physical 
information. 

Our goal in the remaining sections will be to compute the functions T e (a, d) and Gt( a, d) 
for a holographic theory and for a weakly coupled field theory and compare the results. 2 

3 Free Field Theory Results 

We begin by analyzing the case of a free conformal field theory, specifically a massless scalar 
held, in 4+1 dimensions on the background (2.1), following [5]. We will take x c non-conrpact 

2 Since the energy density is simply related to the temperature in a holographic theory, these functions 

carry the same physical information. However, as we shall describe, F e is more natural at weak coupling, while 

Gt is a better diagnostic at strong coupling. 
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to compare with the strongly-coupled results, which are ^-independent. The scalar action is 
in general d + 2 spacetime dimensions is (nb: we will eventually take d = 3) 3 


S = J d d+2 x d v <f>. 


(3.1) 


Expanding (ft in Fourier modes with dk^ = ^ AirS 3 ' 


<Kx,t) = J dk d Mt) e< ts+kcXc ), 


(3.2) 


we obtain a complex harmonic oscillator for each pair {(k, k c ), — (k, A: c )}, with action: 

' k 2 


S = dt a(t) d < | dt(f>k\ - 


- + k 2 

[ a(t ) 2 c 


(3.3) 


These oscillators are completely decoupled, so their quantum states evolve independently. 


3.1 Time-dependent oscillator 

Classically, the mode with momenta ( k , k c ) evolves as 

4> + d '—(f) + ^ 'k 2 -|—2 ^ 4 1 = 0 j (3-4) 

For early and late times, this corresponds to a simple harmonic oscillator, with frequencies 

OJi = sjk 2 + kl UJf = \j +k 2 

respectively. Defining annihilation operators a, b and A, B associated with the oscillators 
at early and late times, we can write the field operator at arbitrary times in two equivalent 
ways 


<£(*) 


a (j>i(t) + bt <%(t ), 
Acj) f (t) + B f <^(t), 


where for g £ {*, /} one has 


—> 



-lUJr! t 
C 1 


t -4 Too. 


(3.5) 


(3.6) 


These are normalized such that <j>i cj>* — (ft* <p t = i. 

3 We could also consider a massive scalar field in d + 1 dimensions, instead of compactifying a massless 
field from d + 2 dimensions down a Scherk-Schwarz circle. However, we choose to perform the computation in 
closer analogy with the holographic set-up for ease of comparison. The massive field answer can be recovered 
by focusing on a single mode in the compact direction, i.e., picking an appropriate value for k c below. 
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(a) (b) 

Fig. 1 : Free field results for: (a) the energy production captured by Afcpp) and (b) the temperature change 
Gt(o, d), as function of the amplitude for a range of o as indicated. This should be contrasted with the strong 
coupling result displayed in Fig. 2. 


As 4>f and form a basis for solutions to (3.4), we have (nb: |ct| 2 — \/3\ 2 = 1) 

= outfit) + P4>*f{t) . (3.7) 

Correspondingly, we have the Bogoliubov transform 

a = a*A - /3* B+, A = aa + 0*b t 

, , (3.8) 

b = a*B-/TA t , B = ab + /3*a t 


If we start with an in-mode in the thermal state 


= —p-P‘ tH = 1 - / 8 T w i (a t a+&t&+l) 


z z 

we can determine the final energy in this mode relative to the vacuum energy. One finds 
E - E vac = TV (u f (At A + Bt B 


(3.9) 


= 2 Uf 


r,~P T Wi 


l — e -/3 T Ui 


+ 


Pt 

2 1 + e _ftrWi 

1 _ g-/3 TWi 


(3.10) 


Here, the first term corresponds to the adiabatic result. We can interpret it as having the 
same occupation probabilities as the initial state, but with modified energies. The second 
term, which vanishes in the limit of slow expansion where the Bogoliubov coefficient j3 goes 
to zero, captures the effects of particle creation. 


3.2 Results 

We can now add the energies for all the modes of the scalar field to obtain an expression for 
the final energy density. The density of modes in k space with an IR regulator L is {^) d+l , 
while the total volume after expansion is L d+l a d . The total energy density is then the integral 
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over all (k,k c ) of (3.10). Normalizing the result correctly by the mode density and the final 
proper volume, we find the final answer (independent of the IR regulator) to be: 1 

e E (a,v,T) = ^ J dk d — _ (3.11) 

For a d = 3 free scalar the initial energy density (equivalently, the result for a = 0) is 

eo(T 0 ) = C^lcTl , Cj, = (3.12) 

Thence from (2.4) and (2.3), we find 

F - ( “' t ’ > = c^/ dk3 “' l/3| -- 1 £^T' <3 ' 13) 

To evaluate F e , we calculate the Bogoliubov coefficient /3(v,k,k c ,a) numerically by finding 
solutions of (3.4). We pick initial conditions 4>i = e~ luJit for early times and read off /3 from 
the decomposition (3.7) in the late time behavior. Our results for the function F e for various 
choices of o are shown in Fig. la. 

Similarly, we can calculate the function Gt{o-, d) defined in (2.5). Note that in the free 
field case, the final state is not thermal. However, assuming some very weak interactions that 
thermalize the system on a time scale much longer than the expansion time, we can define the 
final temperature in terms of the final energy density using the equilibrium relation (3.12). 
Our results for the function Gt(o., d) are shown in Fig. lb. 

4 Holographic Strong Coupling Results 

We now derive the corresponding results for the case of a strongly-coupled holographic field 
theory. For this purpose, we seek asymptotically AdSg solutions to Einstein’s equations with 
negative cosmological constant, with boundary geometry given to be the FRW cosmology 
( 2 - 1). 4 5 

The bulk metric, written in ingoing Eddington-Finkelstein or Bondi-Sachs coordinates, 
is (CvdS = 1): 

ds 2 = -2 He 2 * dt 2 + 2e 2 * dt dr + Y?[e B dx 2 + e~ 3B dx 2 c ) (4.1) 

Since we are interested in homogeneous cosmologies, we take the metric functions {A, x, E, B} 
to depend in the holographic radial direction r and the time coordinate t, but to be indepen¬ 
dent of the spatial coordinates {x,x c }. The boundary scale factor behaves as in (2.2). 

In order to solve the Einstein equations we start the system at an initial thermal state. 
The final state of the system will be also be thermal (even starting from vacuum), and we can 
read off the final energy density, entropy, and temperature of the final state from the solution 
at late times. We can then calculate the quantities F e (a, d) and Gy (a, d) for comparison with 
the free field results. Details of the numerical calculation are given in the Appendix. 

4 Since we integrate over k and —k independently, we have to divide the result of (3.10) by half. 

5 For the strong coupling calculations, we can easily generalize to allow spatial curvature, replacing dx 2 —> 
ds% K , for spatial slices with spatial curvature k = 0, ±1. However, we will mostly discuss the k = 0. 
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4.1 Equilibrium and Adiabatic Physics 


To compute F e and Gt, we need the equilibrium results for the holographic theory. In the 
deconfined phase, the CFT 4 + i compactified on a spatial circle is dual to the Schwarzschild- 
AdSg solution compactified on S, from which we can read off the thermodynamic data: 


e = Ap = AC H T 5 , 


Ch = C e ff 




(4.2) 


The physical quantities on the FRW 4 universe are obtained simply by multiplying these 
expressions by l c . For k ^ 0 we can obtain the equilibrium thermodynamics numerically. 
Furthermore, for an adiabatic expansion from an initial scale factor a* = 1 to a final scale 
factor df = a, keeping £ c fixed, the total entropy should remain fixed. Hence from (4.2) we 
obtain 



CH 


(Too)f 
(Too )i 



(4.3) 


as the proper volume of the FRW 4 scales as a(t) 3 . This dilution factor in energy, which is 
determined by the underlying scale invariance of CFT 44 . 1 , is part way between that for matter 
and that for radiation. This is sensible, since we should have a factor of cH 3 for the volume 
dilution, but the red shift effect operates only in 3 of the 4 spatial directions. 


4.2 Results 

Using the equilibrium (4.2) and adiabatic (4.3) results, we have calculated F e (a, d) and 
Gy (a, 0) numerically for various values of d; see Fig. 2. From Fig. 2a, we note that as 
the amplitude of expansion increases, the function F e (a, d) describing the final energy density 
approaches the same approximately constant curve for each value of the initial temperature 
(or d), in contrast to the free field theory results. This suggests that at strong coupling, 
the field theory quickly forgets about the initial state before the expansion, with the energy 
density in the final state sensitive only to the expansion rate. 

For G(a, d), we also find a convergence of the results for various values of d, but this 
time for sufficiently small a., cf., Fig. 2b. As we will see now, this may be understood as a 
consequence of the fact that in this regime, the strongly-coupled field theory dynamics is well- 
described by relativistic hydrodynamics. Generally, we expect this to be true in the case that 
derivatives are small compared to the scale set by the temperature. In our homogeneous setup, 
physical quantities vary only in time, and the scale of this time dependence is characterized 
by the Hubble parameter H = a/a. Thus, we expect hydrodynamics to be valid when 
H -C T. In this fluid-dynamical regime, the field theory stress tensor is determined by the 


-7 - 




(a) (b) 

Fig. 2: Holographic results at strong coupling as function of the amplitude for a range of o. Contrast against 
Fig. 1 for free field analog. 

(a) The energy production captured by F e (a, o). This quickly becomes insensitive to the initial state as amplitude 
is increased. 

(b) The temperature change Gt(o,d), displaying good agreement with the hydrodynamics prediction (dashed 
line) for small enough c. 


local temperature T(x) and velocity u /x (x), see [6]: 


T > _ /O rpAi 

inv — J- 


T (g^ v + 5 u,j, u v ) 


5 

2n 


P^a.Pv/3 ~ | P^vPa/3 I V“f/ 


+ 0(V 2 


(4.4) 


where P^ u = g fiu + u^u v defines a projection to the spatial directions. The dynamics of T(x) 
and u ,J, (x) are then determined by the conservation relation \7^T^ U = 0. In our case, the 
spatial velocity vanishes, so = 5^. The evolution equation is simply 


■ 3 a 3 

T+-T-= - 

4 a 32vr 



(4.5) 


Starting from initial temperature To with scale factor (2.2), we find the solution Gt(i i, t>) = 
T(t = oo) - 


G T (a, o) 


1 

7tt 


7 

fl 4 


+ 7 a 4 —7a—1 


a 4 (a - 1) 


(4.6) 


The salient feature is that Gt is independent of 0 in the hydrodynamic regime. From 
the numerical results, Fig. 2b, we see very good agreement with this prediction for small n, 





















just as expected from the criterion H <C T. This serves as a strong check of the numerical 
methods. Note that for large enough temperatures the holographic results are well-described 
by hydrodynamics for an increasing range of amplitudes. At lower temperatures, and for 
sufficiently large amplitudes, the non-linear evolution deviates strongly from hydrodynamics. 
It is curious that such deviations are always positive, i.e., lead to more energy and entropy 
production than in hydrodynamics. 

5 Discussion 

In this paper, we have presented results for the evolution of homogeneous states of strongly- 
coupled confining gauge theories in FRW cosmologies. While we have focused on the flat 
k = 0 case, explicit results have also been obtained in the case of positive (k = 1) and 
negative (k = — 1) spatial curvature. We leave this exploration to future study. 

We have found significant qualitative differences from the case of free field theory. First, 
for the range of temperatures considered, the final energy density quickly becomes insensi¬ 
tive to the initial temperature as the amplitude of the expansion becomes large (for fixed 
expansion rate). Second, the small amplitude results are well described by a hydrodynamic 
approximation in which the difference Tf — Tf | adi . lt - )a)ic is independent of the initial tempera¬ 
ture. 

Using similar techniques, it would be straightforward to consider other cosmological 
spacetimes, or to compute other observables such as correlation functions and entanglement 
entropies. A more challenging line of research would be to study the time-dependent decon¬ 
finement phase transition, caused by the expansion. 
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A Appendix: Solution Method 

Equations Of Motion: In choosing the Bondi-Sachs form of the metric (4.1), we are able 
to use the characteristic formulation of Einstein equations. The numerical scheme we use is 
described in great detail in [7], and some of our choices are more similar to the ones made in 
[8]. Here we briefly comment on some issues specific to the models discussed here. We refer 
the interested reader to [7, 8] for a more complete discussion. 
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The Einstein equations, in the nested form used for their solution, are as follows: 

- E B' 2 - 2 S' V + E" = 0 
4 

where prime indicates a radial (r) derivative. We choose to gauge fix the determinant of 
the spatial metric E and solve for the field x at each time step. The gauge fixing is not 
complete: we leave enough freedom to fix the coordinate location of the apparent horizon, for 
convenience. The gauge freedom is implemented in terms of a gauge parameter A (t) which 
we treat as a dynamical variable. 

Next, we solve for the field d+E, where d+ = dt + A d r : 

d + S' + ^±!^-5 e 2xE = 0 (A.l) 


Since the field E is gauge fixed, leaving only the parameter it A in its stead, d + E should be 
thought of as a proxy for the field A and the time derivative of the gauge parameter A, which 
we solve for later after gathering more information (dot denotes time derivative). 

The next equation is for the derivative of the dynamical field B: 


d+B ' + TkT + :l±E£ = „ 


(A.2) 


Once this is solved, we can find A by the requirement that the apparent horizon (defined as 
the locus of the outermost zero of d + E), stays at fixed radial location 


A + - e~ 2x d+B 2 = 0 
6 


When A is obtained, the expression for d+E is sufficient to find the field A, and that 
information in turn can be used to convert knowledge of d+B into an expression for the time 
derivative B. 

The above nested scheme can be used for constrained evolution: at each time step we 
are given the values of the propagating fields {B, A}. We use the above process to determine 
the constrained fields {%, A} as well as the time derivatives of the dynamical fields, {S,A}, 
which are then used to propagate them to the next time step. 

Asymptotic Expansion and Observables: Denote the metric elements in our ansatz 
collectively as gij(t,r). For each such metric element, we can write an asymptotic near¬ 
boundary expansion, by transforming the familiar Fefferman -Graham expansion to incoming 
coordinates. Since in our case (AdS^+i with d odd) there are no logarithmic terms, the 
expansion is simple: 6 
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Fig. 3: Convergence of the total entropy production with the grid size, used to discretize the radial direction. 
Plotted is the (base 10) logarithmic change of total entropy as we refine the grid, for fixed parameters. This 
demonstrates the expected exponential convergence. The plots in the paper were generated with 30 grid points, 
where the error corresponds roughly to machine precision. 


As is well known, all the expansion coefficients up to, but not including are 

determined in terms of the boundary metric g^\t). Since the coefficients are known, we 
choose to shift the fields by these expressions, and solve for the shifted dynamical fields. Such 
choice of variables is described in [7], but in our case the shift includes all the expansion 
coefficients determined by the boundary metric. As a practical matter, since the expressions 
for the expansion coefficients are increasingly complex as the dimensionality increases, the 
resulting equations and boundary conditions are extremely long and uninformative, and we 
spare the reader the precise details. 

To extract the energy-momentum tensor from the asymptotic expansion we use the stan¬ 
dard expressions for the counter-terms and renormalized energy momentum tensor, given e.g., 
in [9]. 

Equilibrium thermodynamics: To illustrate the above, consider the derivation of the 
equilibrium thermodynamic data quoted in (4.2). The Schwarzschild-AdS6 solution with a 
flat boundary metric (k = 0 ) is given by: 


ds 2 


T2 dt2 + dz2 + 9(z) i ds i,o + dw2 )) 

(1 ~Uz 5 ) 2 


(1 + I 


g(z) = (1 + 


(A.3) 


The CFT .5 stress tensor is given by (T^ u ) = 16 Jg . rffi, where rjfj is the z 5 term in the 
asymptotic expansion of the metric. Setting c e g = 16 Jq n as the central charge of the CFT 

6 A similar expansion in situations with d even would lead to logarithmic terms which would appear at 
order ^ log(r). 
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and integrating over the compact S, 1 ,, we have 


(Too) — 4 c e ff n , (Tjj) — c e s i-i 6{j . (A.4) 

The entropy density s computed from the horizon area using the Bekenstein-Hawking formula 

1 4 

S = 4Gv , gives s = 47 rceff/rs (nb: horizon is located at the zero locus f(z + ) = 0). The 
temperature of the black hole (say using dE = T dS ) is obtained to be T = /j,s. 

Note that one of the advantages of working in AdSg is that we do not have to worry 
about logarithmic terms in the asymptotic expansion, which plague odd-dimensional AdS 
spacetimes. 

Numerical Choices: The solution method consists of constrained evolution, which involves 
the iterative solution of several linear ordinary differential equations at each time step. We 
discretize those differential equations using pseudo-spectral collocation methods. To evolve 
the system in time we use a fourth order Runge-Kutta method with an adaptive step size. 
In order to avoid numerical instabilities we use filtering at each time step; we use both filters 
based on Chebyshev interpolation, or filters based on fast Fourier transform. 

Since the expressions involved in the equations are quite long, we need to employ some 
special tricks to minimize round-off errors. To this end we use compensated summation to 
evaluate the sums involved in our equations, and iterative refinement in solving the linear 
equations. Both those steps utilize extra precision in intermediate steps of the calculation. 

We demonstrate the convergence of our solution in Fig. 3. Similar tests were performed 
for other numerical parameters, such as the tolerance involved in determining the temporal 
step size. 
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